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1.  INTRODUCTION 


In  recent  work  (Powell,  Walbert,  and  Zielinski  1993),  we  have  developed  a  model  and 
formalism  for  studying  current  and  heat  transport  in  solid-armature  railguns.  The  model 
is  two-dimensional  and  transient  and  can  be  used  to  determine  the  electromagnetic  fields 
and  forces,  the  current  density,  and  the  temperature  in  both  the  rails  and  the  armature  as 
a  function  of  space  and  time.  Such  an  analysis  is  important  in  order  to  guide  the  design  of 
projectiles,  in  order  to  determine  the  limits  under  which  these  armatures  can  operate,  and 
in  order  to  study  the  dynamics  of  the  armature  and  railgun. 

In  our  previous  work,  we  limited  our  study  to  the  case  of  a  small-bore,  simple,  single¬ 
turn  railgun  in  which  the  current  was  carried  from  rail  to  rail  via  a  U-shaped  or  slotted 
armature.  In  the  present  calculations,  we  will  extend  the  model  to  treat  the  railgun  and 
armature  recently  developed  in  the  Cannon-Caliber  Electromagnetic  Gun  (CCEMG)  Pro¬ 
gram  (Kitzmiller  et  al.  1994).  The  railgun  has  a  rectangular  bore  cross  section  and  is  series 
augmented;  the  armature,  which  is  discussed  in  more  detail  subsequently,  is  considerably 
more  complicated  than  that  which  we  have  previously  analyzed. 

A  schematic  diagram  of  a  series- augmented  railgun,  indicating  the  general  principles  on 
which  it  operates,  is  shown  in  Figure  1.  Both  sets  of  rails  actually  extend  farther  to  the  left 
and  right  of  the  armature  than  is  shown  in  the  figure,  and  are  connected  so  as  to  produce 
equal  currents  in  the  directions  indicated  by  the  arrows.  This  configuration  will  produce 
an  induction  field  in  the  bore  to  the  left  of  the  armature  which  is,  in  the  two-dimensional 
case,  exactly  twice  as  large  as  that  on  the  right.  As  a  result,  the  total  Lorentz  force  which 
accelerates  the  projectile  is  three  times  that  which  occurs  in  a  simple  railgun  for  any  given 
value  of  the  total  current.  The  subject  of  concern  here  is  to  determine  the  distribution  of 
current,  the  resulting  body  forces,  and  the  generation  and  transport  of  heat  everywhere  in 
both  the  rails  and  armature. 

Similar  types  of  calculations  for  simple  nonaugmented  railguns  and  simple  armatures 
have  been  undertaken  by  various  investigators  in  the  past.  In  early  work.  Young  and 
Hughes  (1982)  investigated  the  current-diffusion  problem  analytically  in  the  steady  state  for 
some  rather  simple  types  of  armature  geometry.  Their  work  provided  considerable  insight 
into  how  the  current  could  be  expected  to  vary  and  the  practical  problems  that  would  be  en¬ 
countered  as  the  projectile  velocity  increased.  Later,  their  work  was  extended  to  a  transient 
analysis,  and  some  limiting-case,  exact  solutions  were  obtained  (Hughes  and  Young  1986). 
Similar  transient,  exact  calculations  have  been  undertaken  by  Nearing  and  Huerta  (1989). 
Probably  the  most  complete  treatment  of  the  problem  is  the  work  performed  by  Long  (1986). 
His  early  analysis  was  somewhat  similar  to  that  of  Hughes  and  Young,  and  succeeded  in  pre- 
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Figure  1.  Schematic  Diagram  of  a  Series- Augmented,  Solid-Armature  Railgun. 


dieting  the  current  distribution  analytically  for  simple  armatures  and  rather  low  armature 
velocities  in  the  steady  state.  In  later  graduate- thesis  work,  Long  (1987)  developed  a  finite- 
element  program  with  which  he  successfully  studied  more  complicated  armatures  at  higher 
velocities.  In  addition,  he  included  the  effects  of  Joule  heating,  temperature-dependent  elec¬ 
trical  and  thermodynamic  properties,  and  heat  transport,  all  of  which  had  been  neglected  in 
the  earlier  analytic  calculations.  In  more  recent  work,  Putley  (1989)  adapted  a,  previously 
developed  code  to  solve  the  rail-armature  interaction  problem,  verified  some  of  Long’s  earlier 
calculations  that  were  done  analytically,  and  extended  calculations  to  somewhat  higher  ve¬ 
locities.  All  of  these  calculations  have  been  two-dimensional  so  that  the  rails  were  assumed 
to  be  infinitely  high,  with  no  variation  in  physical  quantities  along  the  direction  of  the  rail 
height. 

In  addition  to  the  two-dimensional  work,  there  are  important  efforts  underway  to  perform 
diffusion  calculations  in  three  dimensions.  Recently,  Yun  and  Price  (1994)  have  performed 
a  three-dimensional  analysis  on  the  CCEMG  armature.  The  model  which  they  employed, 
however,  does  not  account  for  the  motion  of  the  launch  package  (i.e.,  v  —  0).  Early  three- 
dimensional  work  which  does  account  for  conductor  motion  has  been  described  in  several 
references  (Rodger,  Leonard,  Eastman,  and  Atkinson  1989;  Rodger,  Leonard,  and  East¬ 
man  1991;  Rodger  and  Leonard  1991),  and  additional  efforts  (Hsieh  1994)  are  currently 


underway  at  the  Institute  for  Advanced  Technology  to  develop  a  three-dimensional  finite- 
element  code.  Three-dimensional  calculations,  with  appropriate  models  for  physics  at  the 
rail-armature  interface,  should  ultimately  provide  a  realistic  indication  of  the  current  and 
temperature  distributions  in  railguns.  These  models  can  be  used  to  avoid  many  assumptions 
of  very  limited  applicability  associated  with  the  infinite  rail-height  calculations.  However, 
the  work  is  considerably  more  complicated  than  in  two  dimensions  and  is  rather  preliminary 
at  this  point. 

The  calculations  to  be  undertaken  here  are  most  similar  to  those  carried  out  by  Long. 
We  have  developed  a  two-dimensional,  finite-difference  program  which  solves  the  governing 
equations  for  the  fully  transient  case.  Our  motivation  has  been  to  provide  information 
which  could  be  used  in  the  design  of  projectiles,  as  well  as  to  provide  information  which 
could  eventually  be  used  as  input  in  stress-analysis  codes  which  have  been  developed  to 
investigate  in-bore  dynamics.  As  indicated  previously,  the  specific  calculations  described 
here  will  be  for  the  armature  recently  developed  as  an  integral  part  of  the  CCEMG  Program. 
The  organization  of  the  report  is  as  follows.  In  section  2,  the  model  is  described  and  the 
mathematical  formalism  developed.  In  section  3,  we  describe  the  results  of  some  calculations 
for  this  specific  armature  geometry.  Finally,  section  4  contains  our  summary  and  conclusions. 


2.  MODEL  AND  FORMALISM 

The  basic  equations  that  must  be  solved  numerically  consist  of  the  appropriate  Maxwell 
equations  as  well  as  the  energy  transport  equation  for  the  rails  and  projectile.  It  is  convenient 
to  solve  these  equations  in  a  frame  of  reference  in  which  the  armature  is  fixed  with  its  trailing 
edge  located  at  x  =  0  (see  Figure  1),  and  the  rails  move  in  the  negative  x  direction.  We  will 
restrict  ourselves  to  situations  in  which  the  armature  is  symmetric  about  a  plane  centered 
between  the  rails.  This  plane  is  defined  by  y  =  0.  The  rails  and  projectile  are  infinitely 
extended  in  the  z  direction.  The  governing  equations  to  be  developed  are  solved  in  the 
armature  and  the  inner  rails;  the  details  concerning  the  outer  rails  are  not  considered  except 
in  the  sense  that  they  affect  the  boundary  conditions  on  the  magnetic  induction  as  discussed 
previously. 

Let  J,  B,  and  E  be  the  current  density,  the  magnetic  induction  field,  and  the  electric 
field  intensity.  Let  y,  ac,  a,  and  C  be  the  magnetic  permeability,  the  thermal  conductivity, 
the  electrical  conductivity,  and  the  specific  heat.  Finally,  let  T,  p,  and  v  be  the  temperature, 
density,  and  material  velocity  at  a  given  point  in  the  rail-armature  system.  The  governing 
Maxwell  equations  can  then  be  written 
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V  X  5  =  fiJ, 


(1) 


and 


V  X  E  = 


dB 
■  dt  ’ 


(2) 


J  =  (j(E  +  uxE).  (3) 

As  is  customary  in  calculations  of  this  type,  we  have  neglected  the  displacement  current  in 
equation  (1).  This  approximation  can  be  shown  to  be  valid  provided  the  armature  velocity 
is  small  compared  to  the  light  speed  c,  and  provided  the  time  scale  r  of  the  problem  satisfies 
the  condition  r  >>  Both  these  conditions  are  very  easily  satisfied  for  rail-launcher 

problems.  Equation  (3)  represents  Ohm’s  law  in  a  moving  conductor. 

We  now  observe  that  for  the  infinite  rail-height  geometry  discussed  previously,  we  must 
have  that  B  =  Ba^  and  that  E  and  both  vanish.  Furthermore,  there  can  be  no  z 
dependence  in  any  of  the  physical  quantities.  If  we  now  make  use  of  these  observations, 
we  find  that  equations  (l)-(3)  can  be  uncoupled  to  produce  a  single,  second-order,  partial 
differential  equation  which  predicts  the  convection  and  diffusion  of  the  magnetic  induction 
field.  The  equation  can  be  written 

dB  _d^B  d^B  (4) 

~  ^ dx  adTdxdx  adTdydy' 

In  obtaining  equation  (4),  we  have  assumed  that  a  depends  on  position  and  time  only  through 
the  temperature  T, 

A  similar  assumption  for  k  and  C  allows  us  to  write  the  heat-transport  equation  as 


dT  dT  _  1 

dt  ox  fl^cr 


dx)  \dy  ) 


dK  fdTV  dK  (dT\\  d^T  ^  d^T 

i^J 


Equation  (5)  can  be  obtained  from  standard  energy-conservation  principles.  The  meaning 
of  the  various  terms  which  appear  in  the  equation  are  rather  obvious  except  perhaps  for  the 
term  enclosed  in  brackets  on  the  right-hand  side.  If,  however,  we  use  equation  (3),  we  can 
identify  this  term  as  being  equal  to  <7^/cr;  consequently,  it  accounts  for  resistive  heating  in 
the  armature  and  rails. 
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Equations  (4)  and  (5)  are  the  basic  diflterential  equations  which  must  be  solved  subject  to 
some  set  of  initial  and  boundary  conditions.  It  is  evident  upon  examination  of  the  equations 
that,  provided  the  initial  conditions  satisfy  the  conditions  B{y)  =  B{—y)  and  T{y)  =  T{—y), 
symmetry  about  y  =  0  will  persist  for  all  time.  We  will  restrict  ourselves  to  this  case  and 
only  solve  the  equations  for  y  >  0.  We  will  assume  that  there  is  no  transfer  of  heat  from 
the  railgun  or  armature  to  the  atmosphere  surrounding  them.  Consequently,  we  take  for  the 
boundary  conditions  on  T  at  each  conductor  surface  the  relation 


n  •  vr  =  0,  (6) 

where  h  is  the  unit  normal  to  the  surface  (pointing  into  the  conductor).  Such  a  condition 
also  applies  along  the  centerline  y  =  0  from  symmetry. 

Boundary  conditions  on  the  magnetic  induction  can  be  easily  determined  by  arguments 
with  Ampere’s  law.  If  we  denote  by  j  the  current  per  unit  rail  height,  we  must  have  in 
the  inner  bore  to  the  left  of  the  armature  B  =  2fij,  and  to  the  right  B  —  yj.  The  latter 
condition  also  applies  along  the  upper  surface  of  the  inner  rail,  and  along  the  centerline  we 
must  have  from  symmetry  dB/dy  =  0.  At  the  ends  of  the  rail,  we  have  taken  the  condition 
on  the  induction  to  be  dBjdx  =  0.  This  condition  implies  that  the  current  is  parallel  to 
the  rail  surface  at  these  points  and  is  strictly  accurate  only  if  the  ends  are  sufficiently  far 
from  the  armature.  The  validity  of  the  assumption,  as  well  as  an  alternate  approach  for 
handling  these  boundary  conditions,  has  been  discussed  in  previous  work  (Powell,  Walbert, 
and  Zielinski  1993).  Generally,  the  rails  and  armature  consist  of  different  materials  and 
under  such  a  circumstance  dBjdy  and  dTjdy  are  not  continuous  across  the  rail-armature 
interface.  The  appropriate  jump  conditions  were  derived  previously  (Powell,  Walbert,  and 
Zielinski  1993)  and  can  be  written 


a  dy 


=  0, 


interface 


(7) 


where  the  brackets  (here  only)  denote  the  change  in  the  enclosed  quantity  as  the  interface 
is  crossed.  A  similar  relation  for  the  temperature  holds  if  B  is  replaced  by  T  and  1/a  by  k. 


We  finally  consider  the  calculation  of  the  induction  field  in  regions  of  the  armature  which 
are  bounded  on  all  sides  by  conducting  surfaces.  Such  a  region  is  indicated,  for  example, 
between  the  front  and  back  contacts  of  the  armature  in  Figure  1.  Since  the  conductivity  in 
this  region  is  zero,  it  follows  from  equations  (l)-(3)  that  B  here  is  independent  of  position. 
We  denote  the  constant  value  by  Bh,  the  area  of  the  region  by  Ah,  and  integrate  equation  (2) 
over  the  area.  Application  of  Stokes’  theorem  then  yields  the  relation 
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(8) 


Ah 


dB 

dt 


=  -/- 


■dl. 


where  the  integration  denotes  a  line  integral  in  the  counterclockwise  direction  around  the 
periphery  of  the  defining  area.  In  obtaining  this  result,  we  have  used  equation  (3)  as  well  as 
the  fact  that  the  tangential  component  of  E  must  be  continuous  along  the  line  of  integration. 
The  constraint  imposed  by  equation  (8)  is  sufficient  to  determine  the  time  evolution  of  the 
spatially  constant  induction  Bh- 

Equations  (4),  (5),  and  (8)  with  the  indicated  boundary  conditions  are  sufficient  for  de¬ 
termining  all  electrodynamic  properties  as  well  as  the  temperature  as  a  function  of  both 
space  and  time  within  the  rails  and  the  armature.  We  have  previously  developed  an  im¬ 
plicit,  finite-difference  code  for  solving  the  relevant  partial  differential  equations.  This  code 
has  now  been  extended  to  treat  augmented  railguns  and  reasonably  complicated  armature 
geometry.  The  details  of  the  numerics  have  been  discussed  elsewhere  (Powell,  Walbert,  and 
Zielinski  1993)  and  are  not  repeated  here. 


3.  CALCULATIONS 

We  now  discuss  two  calculations  that  we  have  done  to  analyze  current  and  heat  transport 
in  the  armature  developed  by  the  University  of  Texas  Center  for  Electromechanics  (UT- 
CEM)  in  the  CCEMG  Program.  The  planar  cross  section  of  the  armature  appears  as  shown 
schematically  in  Figure  1.  This  armature  is  designed  to  reach  a  launch  velocity  of  1.85  km/s 
via  application  of  a  half  cycle  of  a  sinusoidal  current  pulse  applied  to  a  series-augmented 
railgun.  The  pulse  has  a  duration  of  2  ms  with  a  peak  current  of  780  kA.  In  the  first  case 
undertaken  here,  the  armature  is  accelerated  to  1.85  km/sj  in  the  second,  the  armature  is 
held  fixed  so  that  we  can  examine  directly  the  effects  of  velocity  upon  the  resulting  magnetic 
induction  and  temperature  profiles. 

As  is  well  known,  two-dimensional  models  directly  applied  severely  overestimate  the 
Lorentz  forces  and,  consequently,  the  acceleration  imparted  to  the  projectile  even  in  sim¬ 
ple  railguns.  The  effect  is  considerably  exacerbated  in  augmented  guns.  There  have  been 
developed  some  rather  ad  hoc  approximations  (Batteh  1984)  for  reducing  the  magnitude  of 
these  forces,  and  these  approximations  work  reasonably  well  for  simple  railguns  in  the  steady 
state.  They  cannot,  however,  be  expected  to  work  in  the  transient  case.  In  an  effort  not  to 
overestimate  both  the  Lorentz  force  and  the  Ohmic  heating  at  early  times,  we  have  scaled 
the  current  per  unit  rail  height  down  appropriately  to  obtain  the  correct  acceleration  for  a 
truly  two-dimensional,  infinite  rail-height  launcher.  More  specifically,  we  assumed  a  current 
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profile  identical  in  form  to  that  in  the  actual  design,  namely, 

j  =  josin{7dlto),  (9) 

and  then  calculated  the  amplitude  jo  so  that  the  two-dimensional  armature  reached  a  final 
velocity  of  1.85  km/s  in  a  time  U  of  2  ms.  The  value  obtained  for  jo  was  16.1  MA/m. 

Numerical  calculations  were  then  carried  out  for  an  aluminum  armature  and  copper  rails. 
Temperature-dependent  material  properties  for  these  metals  have  been  given  in  our  earlier 
report  (Powell,  Walbert,  and  Zielinski  1993).  We  now'  discuss  some  of  the  more  significant 
results  of  the  calculations. 

Shown  in  Figure  2  are  lines  of  constant  magnetic  induction  a,t  1  m.s  (peak  current)  and  at 
the  termination  of  the  pulse' at  2  ms.  Only  the  top  half  of  the  armature  is  shown,  since  there 
is  symmetry  about  the  y  a.xis,  The  corresponding  armature  velocities  at  these  times  were  925 
and  1,850  m/s.  It  is  possible  to  prove  from  Maxwell’s  equations  that  current  cannot  cross 
lines  of  constant  induction  so  these  lines  also  represent  streamlines  of  the  current  profile.  At 
1  ms,  when  the  total  current  is  maximum,  a  phenomenon  known  as  the  velocity  skin  effect 
is  very  appai'ent  from  the  curvature  of  the  field  lines  tow'a,rd  the  left  at  the  rail-armature 
interface.  In  the  back  contact,  the  current  is  directed  in  the  positive  y  direction  and  very 
much  concentrated  toward  the  back  of  the  contact.  This  concentration  of  current  arises 
at  even  moderate  velocities  (Long  1986)  because  there  is  insufficient  time  for  current  to 
diffuse  into  the  rails  during  the  time  the  armature  is  in  contact  with  a  given  point  on  the 
rail  surface.  For  example,  we  can  crudely  estimate  from  the  well-known  skin-depth  formula, 
5  =  (TrtliicrY^'^,  that  for  a  stationary  armature  current  will  diffuse  into  the  rails  a  distance 
of  approximately  1  cm  during  the  2  m,s  time  of  the  current  pulse.  Consequently,  the  mean 
“diffusion  velocity”  is  only  about  5  m/s,  and  arro.ature  velocities  of  only  a  few  tens  of 
meters  per  second,  new  rail  is  being  encountered  too  fast  for  the  current  to  diffuse  enough  to 
produce  any  appreciable  skin  depth.  Consequently,  the  resistance  along  the  interface  remains 
high,  and  the  current  remains  confined  to  the  corner.  In  the  front  contact,  the  current  is  in 
the  opposite  direction,  and  here  the  velocity  tends  to  enhance  the  transport  of  current  across 
the  contsxt  at  points  near  the  interface.  At  points  far  from  the  interface,  hovv'ever,  the  current 
is  still  largely  confined  to  the  front  and  back  sides  of  the  armature  at  this  fairly  early  time. 
At  2  ms,  w'hen  the  total  current  has  vanished,  considerably  more  diffusion  into  the  armature 
has  occurred.  It  should  be  noted  that  there  is  significant  magnetic  energy  in  the  armature 
because  of  eddy  currents  even  at  the  termination  of  the  pulse.  This  energy  results  primarily 
from  the  fairly  rapid  rise  time  and  decay  of  the  sine  wave.  It  should  also  be  observed  that 
the  hollow'  region  in  the  armature  between  the  contacts  is  effective  in  enhancing  the  rate  of 
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diffusion  across  the  armature  since  this  region  has  negligible  conductivity.  The  induction 
field  in  this  region  is  not  shown  in  the  figure. 

Shown  in  Figure  3  are  isotherms  at  2  ms,  the  termination  of  the  current  pulse,  when 
the  final  armature  velocity  was  1,850  m/s,  Shown  for  comparison  are  results  of  an  identical 
calculation  except  that  the  armature  was  arbitrarily  held  at  rest  (u  =  0).  No  phase  change 
was  taken  into  account  in  the  model  even  though  at  some  few  points  the  melting  temperature 
of  aluminum  (~  930  K)  was  exceeded.  Here  again  the  velocity  skin  effect  is  quite  evident 
from  the  substantially  greater  heating  that  occurs  at  the  left-hand  edge  of  the  rail-armature 
interface  (i.e.,  the  rear  contact)  for  the  calculation  in  which  the  velocity  was  accounted  for. 
Clearly,  there  is  significantly  greater  heating  on  the  left  side  of  the  armature  than  on  the  right 
because  of  the  much  larger  eddy  current,  particularly  at  early  times.  Toward  the  interior  of 
the  armature,  the  current  is  sufficiently  small  and  of  sufficiently  short  duration  that  negligible 
heating  occurs.  It  should  also  be  observed  that  there  is  enhanced  heating  at  sharp  corners  in 
the  armature  where  the  angle  measured  in  the  conductor  is  greater  that  180  degrees.  This 
effect  has  been  noted  in  our  previous  work  (Powell,  Walbert,  and  Zielinski  1993). 

Further  demonstration  of  the  importance  of  the  velocity  skin  effect  can  be  seen  in  Figure  4 
in  which  the  force  loading  in  the  armature  is  plotted  at  peak  current  (t=l  ms)  for  both  u  =  0 
and  u=925  m/s.  The  force  density  indicated  is  simply  the  magnitude  of  JxB.  At  local  points 
this  force  has  both  x  and  y  components  but,  when  averaged  over  the  armature,  produces  a 
net  result  which  has  only  an  x  component  and  accelerates  the  projectile  down  the  barrel  of 
the  gun.  As  is  evident  in  the  figure,  the  large  current  concentration  at  the  left-most  corner 
of  the  rail-armature  interface  when  v  ^0  produces  significantly  greater  forces  in  that  region 
than  occur  in  the  u  =  0  case.  Similarly,  the  forces  in  the  right  contact  near  the  interface  are 
larger  when  u  /  0  because  of  the  transport  of  current  across  the  contact  near  the  interface. 
At  points  far  removed  from  the  interface,  the  force  distribution  looks  much  the  same  in  the 
two  cases  as  does  the  temperature  distribution  for  the  two  cases  discussed  in  Figure  3. 

To  illustrate  more  qualitatively  some  of  the  previously  mentioned  observations  we  have 
plotted  in  Figure  5  the  total  magnetic  energy  per  unit  rail  height,  Em,  in  the  armature  and, 
in  Figure  6,  the  current  per  unit  rail  height  in  the  front  and  back  contacts.  The  quantities 
are  given  as  functions  of  time  throughout  the  course  of  the  current  pulse.  For  both  plots 
we  present  results  which  account  for  the  velocity  of  the  armature  as  well  as  results  for  the 
calculation  in  which  the  armature  was  held  at  rest.  The  magnetic  energy  per  unit  rail  height 
was  calculated  from  the  expression 


Em  = 


(10) 


9 


Figure  3.  Isotherms  in  the  Armature  at  2  ms 


with  the  integration  being  taken  over  the  planar  cross-sectional  area  of  the  armature.  The 
current  per  unit  rail  height  through  the  contacts  (directed  in  the  y  direction)  can  be  obtained 
by  integrating  equation  (1)  from  the  trailing  to  the  leading  edge  on  the  contact  in  question. 
If  we  normalize  the  resulting  values  by  the  total  current  per  unit  rail  height  j,  we  have 


jp  -  jp/j  -  Bn/jJ-  -  1 

-  Wi  =  2  -  (11) 

where  F  and  R  denote  the  front  and  rear  contacts,  respectively. 

As  i,s  evident  in  Figure  5,  the  peak  magnetic  energy  lags  the  peak  current  (not  shown), 
which  occurs  at  1  ms  by  several  hundred  microseconds.  Furthermore,  there  i,s  a  significant 
amount  of  magnetic  energy  (c^i  45  kJ/rn)  stored  in  the  arm.ature  after  the  current  pulse  has 
ceased.  This  energy  arises  from  eddy  currents  which  i)i  turn  are  caused  by  the  inability  of 
the  armature  to  respond  to  the  current  in  a  quasi-stationary  mannei'.  At  early  times  when 
the  velocity  is  small,  both  curves  in  Figure  5  yield  essentially  the  same  results.  At  later 
times  the  magnetic  energ}''  is  slightly  higher  wdien  v  ^  0  compared  to  the  no-velocity  case. 
Presumably,  the  higher  value  results  from  the  enhanced  transport  of  the  augmentation  field 
across  the  interface  from  right  to  left  when  v  ^  0.  V/e  have  actually  carried  the  calculation 
out  an  additional  250  y.s  with  j  =  0,  i.e.,  to  2,250  /is,  for  the  no-velocity  case  and  found 
that  the  energy  had  decayed  to  about  15  kJ/m  at  that  time.  This  value  represents  about 
5%  of  the  peak  value.  Although  the  amount  of  magnetic  energy  stored  in  the  armature  at 
the  temiinaiion  of  the  current  pulse  is  relatively  small,  the  resultant  defiections  and  stresses 
may  perturb  the  arrna.ture  and  affect  the  discard  of  the  sabot.  This  phenomenon  is  beyond 
the  scope  of  this  study  and  requires  further  investigation. 

In  Figure  6  are  sb.ovtrn  the  normalized  current  per  unit  rail  height  in  the  the  front  and 
back  contacts,  jp  and  as  a  function  of  time.  The  neg.a.tive  signs  are  indicative  of  current 
directed  in  the  negative  y  direction  and,  of  course,  =  1.  At  early  times,  the  magnitude 

of  the  current  through  the  back  contact  is  exactly  twice  that  through  the  front.  This  behavior 
results  from  the  larger  induction  field  at  the  back  and  the  negligible  amount  of  current 
tra.nsport  across  the  armature.  It  is  particularly  noteworthy  that  for  the  ca,se  of  no  velocity 
the  current  distribution  is  fairly  constant  during  the  time  of  the  current  pulse.  However, 
when  the  velocity  is  accounted  for,  the  transport  of  current  from  the  front  toward  the  back 
surface  leads  to  a  m.ucb  more  uniform  distribution.  In  fact,  at  the  very  latest  times  there 
is  nearly  an  equal  amount  of  current  through  the  front  and  back  contacts,  and  both  are 
positively  directed.  Irrespective  of  direction,  the  current  partition  at  the  back  and  front 
contacts  is  roughly  2:1  for  the  v  =  0  case.  For  the  case  with  velocity',  the  back  and  front 
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contact  current  partition  is  more  time-dependent  and  reaches  a  maximum  of  4:1  shortly  after 
peak  current. 

In  Figure  7  are  shown  a  series  of  isotherms  in  the  rail  at  500- /zs  intervals.  For  presentation, 
we  have  transformed  these  results  to  the  reference  frame  in  which  the  rails  are  fixed.  That 
transformation  is  accomplished  via  the  relation  X  =  x  —  Xend-,  where  x^nd  denotes  the  location 
of  the  left-most  end  of  the  rail  at  the  time  in  question.  The  armature,  not  shown  in  this 
figure,  is  located  very  near  the  right-most  end  of  the  rail  (see  Figure  8);  the  part  of  the 
rail  farther  to  the  right  is  excluded  from  the  figure.  There  is  some  heating  at  the  top  of 
the  rail  resulting  from  the  diffusion  of  the  augmenting  field  {B  =  fij  at  the  upper  surface) 
downward.  More  significant,  however,  is  the  heating  of  the  inner-bore  surface  of  the  rail.  The 
larger  temperatures  there  result  both  from  the  larger  induction  field  at  the  surface  (B=2  //;), 
as  well  as  the  large  current  concentration  at  the  interface  corner  of  the  armature  contact. 
Clearly,  the  amount  of  heating  depends  upon  the  past  history  of  the  current  profile,  as  well 
as  the  speed  of  the  armature  as  it  passed  the  point  in  question.  In  particular,  high  induction 
fields  produce  high  currents  and  significant  Ohmic  heating;  high  velocities  produce  large 
relative  current  concentration  at  the  interface  corner,  but  small  dwell  times  in  the  vicinity 
of  a  given  point  in  space.  Thus,  at  late  times,  when  the  velocity  is  high  and  the  current  is 
decreasing  rapidly,  relatively  little  heating  near  the  rail-armature  interface  occurs.  It  should 
also  be  noted  that  there  is  insignificant  heating  in  the  interior  of  the  rail.  This  result  suggests 
that  there  is  incomplete  penetration  of  the  current  across  the  rail  thickness,  a  distance  of 
2  cm. 

We  finally  show  in  Figure  8  an  enlarged  view  of  the  rail  isotherm  at  peak  current.  Only 
the  first  5  mm  above  the  bore  surface  of  the  rail  is  shown.  (It  should  be  noted  that  the 
armature  is  not  drawn  to  scale  in  the  y  direction.)  More  detail  is  evident  than  exists  in 
the  previous  figure,  and  the  profiles  closely  resemble  those  discussed  in  our  previous  work 
(Powell,  Walbert,  and  Zielinski  1993).  It  is  noteworthy  that  the  maximum  temperature  at 
the  rail  surface  is  only  about  700  K,  and  thus  lower  than  both  the  melt  temperature  of 
copper  and  the  temperature  at  “hot  spots”  in  the  armature.  While  this  temperature  rise 
is  of  minor  concern  for  the  single-shot  computations  discussed  here,  CCEMG  is  a  multishot 
system,  and  the  resulting  long-term  temperature  effects  need  further  study. 


4.  SUMMARY  AND  CONCLUSIONS 

We  have  employed  a  two-dimensional  model  to  investigate  current  and  heat  transport  in 
an  armature  developed  by  UT-CEM  in  the  CCEMG  Program.  The  armature  is  accelerated 
to  1.85  km/s  in  a  series-augmented  railgun.  Time-dependent  calculations  were  carried  out 
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Note:  Armature  not  drawn  to  scale  In  the  y  direction. 


for  the  duration  of  the  current  pulse  with  the  armature  motion  accounted  for  and  compared 
with  calculations  in  which  the  armature  was  held  fixed.  General  conclusions  reached  as  a 
result  of  this  analysis  can  be  reiterated  as  follows. 

(1)  The  velocity  skin  effect  is  significant  at  points  near  the  rail-armature  interface  at  even 
moderate  velocities.  At  points  far  removed  from  the  interface,  the  effect  is  small  and  the 
induction  and  temperature  profiles  appear  similar  to  those  for  the  no-velocity  calculation. 

(2)  There  is  significant  thermal  and  force  loading  at  the  left-hand  corner  of  the  rail- 
armature  interface  because  of  the  high  current  density  produced  by  the  velocity  skin  effect. 
In  addition,  there  is  significant  heating  and  force  loading  at  large-angle  corners  within  the 
armature. 

(3)  Near  the  rail-armature  interface,  the  transport  of  the  augmenting  induction  field 
across  the  armature  is  enhanced  by  both  the  velocity  skin  effect  and  by  the  hollow  region 
between  the  contacts.  For  the  case  of  no  velocity,  however,  or  at  points  far  removed  from  the 
interface,  there  is  insufficient  time  during  the  shot  for  the  field  to  diffuse  completely  across 
the  armature. 

(4)  The  fast  rise  time  and  decay  of  the  current  pulse  lead  to  some  magnetic  energy 
remaining  in  the  armature  (and  rails)  at  the  termination  of  the  pulse. 

(5)  Thermal  loading  in  the  rails  is  small  relative  to  that  at  localized  points  in  the  arma¬ 
ture.  However,  a  region  of  the  rail  surface  near  the  breech  reaches  temperatures  of  about 
700  K  for  the  single  shot  considered. 

There  are  obviously  significant  limitations  with  the  use  of  two-dimensional  models  for 
quantitative  predictions  in  real  railguns  with  finite  rail  height.  We  do  believe,  however,  that 
these  calculations  are  useful  for  obtaining  a  qualitative  understanding  of  how  these  devices 
operate  and  for  studying  physical  phenomena  associated  with  them.  Primary  among  these 
phenomena  is  the  physics  which  occurs  at  the  rail-armature  interface.  Because  of  the  high 
relative  velocities  at  that  boundary  as  well  as  the  existence  of  a  sharp  corner  on  the  left-hand 
side  of  the  armature,  we  have  been  unable  to  accurately  resolve  the  temperature  and  field 
profiles  near  that  corner.  The  problem  exists  in  a  region  of  only  a  few  grid  points  near  the 
corner  and  is  probably  not  too  important  for  the  problem  studied  here.  However,  in  future 
calculations  in  which  phase  changes  are  accounted  for,  melting  is  likely  to  always  occur  near 
such  a  sharp  corner,  and  it  is  important  to  obtain  accurate  predictions  near  that  point. 
It  seems  likely  that  alternate  numerical  techniques  as  well  as  superior  models  for  handling 
the  interface  physics  are  likely  to  be  required.  We  are  investigating  these  techniques  at  the 
present  time. 
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ATTN  L  HOLLAND 
PO  BOX  85608 
SAN  DIEGO  CA  92138 

1  ELECTROMAGNETIC  RESEARCH  INC 

ATTN  PETER  MONGEAU 
2  FOX  ROAD 
HUDSON  MA  01749 


2  LAP  RESEARCH  INC 

ATTN  JOHN  P  BARBER 
DAVID  P  BAUER 
2763  CULVER  AVE 
DAYTON  OH  45429-3723 

2  SCIENCE  APPLICATIONS  INTL  CORP 

ATTN  K  A  JAMISON 
G  ROLADER 

1247  B  NORTH  EGLIN  PARKWAY 
SHALIMAR  FL  32579 

1  UDLP 

ATTN  BRAD  GOODELL 
MAILSTOP  M170 
4800  E  RIVER  ROAD 
MINNEAPOLIS  MN  55421-1498 

2  SCIENCE  APPLICATIONS  INTL  CORP 
ATTN  JAD  H  BATTEH 

L  THORNHILL 

1503  JOHNSON  FERRY  RD  SUITE  100 
MARIETTA  GA  30062 

1  WESTINGHOUSE  SCIENCE  &  TECHLGY  CTR 

ATTN  DOUG  FKSE 
1310  BEULAH  ROAD 
PITTSBURGH  PA  15233 

1  AUBURN  UNIVERSITY 

ATTN  EUGENE  CLOTHIAUX 
206  ALLISON  LAB 

AUBURN  UNIVERSITY  AL  36849-5311 

1  TEXAS  TECHNICAL  UNIVERSITY 

ATTN  M  KRISTIANSEN 
DEPT  OF  EE  COMPUTER  SCIENCE 
LUBBOCK  TX  79409-4439 

5  INSTITUTE  FOR  ADVANCED  TECHLGY 

ATTN  H  FAIR 
KT  HSIEH 
D  MCNAB 
P  H  SULLIVAN 
J  V  PARKER 

4030  2  W  BRAKER  LANE 
AUSTIN  TX  78759-5329 
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NO.  OF 

COPIES  ORGANIZATION 

1  UNIVERSITY  OF  MIAMI 

ATTN  MANUEL  A  HUERTA 
DEPT  OF  PHYSICS 
PO  BOX  248046 
CORAL  GABLES  FL  33124 

I  NORTH  CAROLINA  STATE  UNTV 

ATTN  JOHN  GILLIGAN 
DEPT  OF  NUCLEAR  ENGRG 
PO  BOX  7909 
RALEIGH  NC  27695 

I  CTR  FOR  LASER  APP  ENGR  SCI  &  MECH 
ATTN  DENNIS  KEEFER 
UNIV  OF  TENNESSEE 
SPACE  INSTITUTE  MAIL  STOP  14 
B  H  GOETHERT  PARKWAY 
TULLAHOMA  TN  37388-8897 

3  CTR  FOR  ELECTROMECHANICS 

ATTN  SIDDHARTH  B  PRATAP 
JOHN  H  PRICE 
HEEDO  YUN 

UNIV  OF  TEXAS  AT  AUSTIN 
10100  BURNET  ROAD  BLDG  133 
AUSTIN  TX  78748 

1  LORAL  VOUGHT  SYSTEMS 

ATTN  ROBERT  TAYLOR 
PO  BOX  650003  WT  17 
DALLAS  TX  75265-0003 


NO.  OF 

COPIES  ORGANIZATION 

ABERDEEN  PROVING  GROUND 

1  CDR  USATECOM 

ATTN  AMSTE  SI  F 

31  DIR,  USARL 

ATTN:  AMSRL-W 

DR.  C.  H.  MURPHY 
AMSRL-WT-WD 
DR  A.  NIILER 
DR  P.  BERNING 
DR  R.  BOSSOLI 
DR  S.  CORNELISON 
DR  A.  GAUSS 
DR  C.  HOLLANDSWORTH 
DR  C.  HUMMER 
MR  L.  KECSKES 
DR  T.  KOTTKE 
DR  M.  MCNEIR 
DR  J.  POWELL  (5  CPS) 

DR  A.  PRAKASH 
MR  C.  STUMPFEL 
DR  G.  THOMPSON 
AMSRL-WT-WE,  DR  J.  TEMPERLEY 
AMSRL-WT-P,  MR.  A.  HORST 
AMSRL-WT-PA 
MS.  G.  WREN 
MR.  G.  KATULKA 
AMSRL-WT-PB 
DR.  E.  SCHMIDT 
MR.  A.  ZIELINSKI  (5  CPS) 
AMSRL-WT-PD 
DR.  W.  DRYSDALE 
MR.  DO  HOPKINS 
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Intentionally  left  blank. 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your  comments/answers 
to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  ARL  Report  Number  ARL-MR-258  _ _ _ Date  of  Report  August  1995 _ _ 

2.  Date  Report  Received  _ _ _ _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for  which  the  report 

will  be  used.) _ _ _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source  of  ideas,  etc.) 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved,  operating  costs 
avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate  changes  to 
organization,  technical  content,  format,  etc.) _ 


Organization 


CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address  above  and  the 
V  Old  or  Incorrect  address  below. 


Organization 


OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  tape  closed,  and  mail.) 
(DO  NOT  STAPLE) 


